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Abstract 

We study the modification of the nucleon nucleon interaction in a nuclear medium 
in the scalar isoscalar channel, mediated by the exchange of two correlated (a chan- 
nel) or uncorrelated pions. For this purpose we use a standard approach for the 
renormalization of pions in nuclei. The corrections obtained for the NN interaction 
in the medium in this channel are of the order of 20% of the free one in average, and 
the consideration of short range correlations plays an important role in providing 
these moderate changes. Yet, the corrections are sizable enough to suggest further 
studies of the stability and properties of nuclear matter. 
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1 Introduction 



The determination of the binding energy of nuclei starting from realistic NN potentials is one of 
the subjects which has received permanent attention from the early days when the Brueckner- 
Bethe-Goldstone (BBG) equation introduced methods to overcome the strong repulsion of the 
nuclear forces at short distances. 

At present, several many body techniques compete to accurately determine the binding energy 
of nuclear matter starting from the realistic NN potentials. One of them is the correlated basis 
functions (CBF) [1,2], which follows the line of the Hypernetted Chain Approach (FHNC) [3]. 
Another one follows the traditional BBG approach [4], and, although costly numerically, the 
variational Monte Carlo method (MC) has allowed to make, in principle, exact calculations, 
although limited to nuclei with small or medium value of A [5] . Methods like the selfconsistent 
treatment of the nucleon selfenergy have also introduced new advances in the field [6-8]. 

The need for three body forces has also been emphasized and the present status is that it is 
difficult to be quantitative on the strength of these forces, and usually they are parameterized in 
order to adjust the precise value of the binding energy [5,9, 10]. It has also been noted in [11-13] 
that short range correlations play an important role when considering these three body forces. 

A common feature of these approaches is that they start from the realistic nucleon nucleon 
interaction, obtained from fits to NN data and deuteron data. They use hence the free NN 
interaction as input. One of the important ingredients of this interaction is the one pion exchange 
(OPE). However, from detailed studies of the pion nuclear interaction it is well known that the 
pion properties in the nuclear medium are sizably renormalized [14-20]. 

There is also the question of the intermediate range attraction, which is basic in the binding 
of nuclei. Models for this interaction would contain a exchange, uncorrelated two pion exchange 
and omega exchange [21]. In as much as the pion properties are changed in the medium, so 
should the two pion exchange be modified. Medium effects in the two-pion exchange have been 
investigated in early works like in Ref. [22] restricting themselves in this case to a subset of two 
pion exchange diagrams with no A-isobar intermediate states, by including Pauli blocking in the 
intermediate nucleons. 

The medium modification of the two pion exchange got a new impulse after the models of 
the 7T7T interaction in the medium showed large modifications [23,24], later on softened by the 
introduction of chiral constraints [25,26]. The implementation of this medium modified tttt inter- 
action in the correlated two pion exchange AW potential increased appreciably the NN attraction 
in nuclear matter. This was partially reduced by the consideration of the chiral constraints in 
Refs. [27,28]. In these latter references the importance of short range correlations which modify 
the 7r nucleus selfenergy was already discussed. The further use of medium modified vector meson 
masses led to improvements in the nuclear matter saturation curve. 

A new perspective into this problem has been made possible by studies of meson meson 
interaction within chiral unitary approaches [30-34] which allow to improve the description of 
the correlated two pion exchange AW interaction [29], as an alternative to the conventional a- 
meson exchange interaction. This picture of the a exchange was mandatory after extensive studies 
showing that the a is not a genuine resonance, made up of qq but just the manifestation of a 
pion - pion resonance state created by the interaction of the pions, what is called a dynamically 
generated resonance. This shows up naturally within the context of chiral unitary approaches 
which use the input of the chiral Lagrangians for the meson meson interaction and extends chiral 
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perturbation theory to implement exactly unitarity in coupled channels [30-34]. This means the 
a exchange inside a nuclear medium will also be modified as direct consequence of the change of 
the pion properties. 

Our aim in the present paper is to start from this new picture for the a exchange, use also 
the standard approach for the uncorrelated two pion exchange and modify these in the nuclear 
medium to see what changes one finds from these sources. Further improvements come from the 
consideration of short range correlations not only in the pion selfenergy but also in the vertex 
functions appearing in the model. The changes obtained are moderate, thanks to the simultaneous 
consideration of these nuclear short range effects in the calculation. In the absence of these, the 
renormalization of the NN interaction is huge. Yet, even the moderate results obtained are large 
enough to motivate further calculations of the nuclear binding and other properties of matter. 

The paper is organized as follows. In Sec. II we provide those elements of the chiral Lagrangian 
which are relevant for the present calculations and briefly discuss peculiarities of the finite baryon 
density. In Sec. Ill we consider the modification of the one pion exchange NN force in the nuclear 
medium and in Sec. IV we discuss the propagation of two pions in the nuclear matter. Sec. V 
and VI are devoted to the in-medium two pion exchange in the scalar-isoscalar channel, both, 
correlated and uncorrelated. The technical details are relegated to the Appendix. 



2 Effective Lagrangian 

In this section we will briefly specify those elements of the effective chiral Lagrangian in the 
meson-baryon sector which are relevant for the subsequent calculations. 

The effective chiral Lagrangian is written as the sum of a purely mesonic Lagrangian Cm and 
the baryonic Lagrangian Lb 

£*ff = £m + (1) 

Both are organized in a derivative and quark mass expansion. The lowest order mesonic La- 
grangian £2 is given by 

C 2 = il(d^d^U + X U ] + U X ] ) (2) 

and contains the most general low-energy interactions of the pseudo-scalar meson octet. In Eq. (j2J) 
the symbol (...) indicates the trace in flavor space, the Goldstone fields are collected in a unitary 
matrix U, f n ~ 93 MeV is the pseudoscalar decay constant and the leading symmetry-breaking 
term \ is linear in the quark masses. For SU(2) and in the isospin limit % = diag(m^,m^). The 
lowest order baryon octet - meson octet Lagrangian reads 

4 1} = (B(i^D, - m B )B) + ^(B^TSK. B V + ^WtsK*, B]) (3) 

where the brackets [...] and {...} denote commutators and anti-commutators, respectively. The 
covariant derivative of the SU(3) baryon matrix B is defined as 

D^B = d, L B + [T M , B] (4) 

In the absence of external field Eqs. © and (jlj) involve other quantities 

u = VU, u> 1 = iuW^Uu^, = ^d^u + ud^v)) (5) 
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Figure 1 : The wNN vertex (a) , the contact itititNN (b) and pion pole (c) terms. 



The SU (3) axial vector coupling constants are determined by neutron and hyperon /3-decay. One 
finds F ~ 0.51, D ~ 0.76 and the axial coupling constant is qa = F + D ~ 1.27. In the SU{2) 
limit the Lagrangian simplifies to 



6$ = ${%>fD 



D + F - 

fj, m N )ip H ~ — ^7 75^^ 



(6) 



where tp is a two component Dirac field ip = (p,n) T . 

In the pion-nucleon sector the chiral Lagrangian (|SJ) constrains all possible interactions of 
the pion fields with fermions at the lowest chiral order we are working. For instance, using the 
exponential parameterization of the unitary matrix U 



U = exp 



T7T 



(7) 



where r are the Pauli operators and expanding u and we obtain the pion-nucleon couplings, 
including up to three pion fields 



i 



U 24/3 

If we supplement Eq. with the lowest order interaction of pions as provided by Eq. (j2J 



(8) 



48/2 



(9) 



we arrive to the set of Feynman graphs shown in Fig. ^ Here, in addition to the standard 
irNN vertex, Fig. the chiral perturbation theory generates the contact term of Fig. ^> (see 
Appendix A) and the pion pole term of Fig. These last two diagrams are the basic element in 
the description of ttN — > ttttN reaction near the threshold [35-38] . The appearance of the pole 
terms and contact 3irNN interactions at the same order of the chiral expansion is crucial for the 
in-medium calculations where due to the partial cancellations the physical amplitudes become 
independent of the parameterization of U matrix even in the presence of the nuclear background 
in accord with the equivalence theorem [39]. One can see explicitely, that the contact term (b) 
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Fi gure 2: The OPEP diagrams. The first graph is the vacuum contribution. The second and third 
diagrams of the upper and first two diagrams of the lower sets of diagrams correspond to the p — h and 
A — h RPA series. The last two diagrams account for the excitation of 2p — 2h states and the contribution 
of the S'-wave optical potential, respectively. The latter diagrams play a minor role and will not be further 
discussed. 



cancels exactly the off shell part of the pion pole term (c) coming from the (q 2 — m 2 ) part of the 
7T7T — ► tttt vertex [40] . We shall also see that the off shell part of the tttt — > tttt amplitude cancel 
exactly with other terms when we perform the calculation of the NN interaction in the nuclear 
medium. 

3 One-pion exchange at finite density 

We start with the one-pion exchange NN potential (OPEP). The typical diagrams modifying it 
are shown in Fig. [2] where the propagation of exchanged pions is distorted by interactions with 
nucleons forming the Fermi sea. The fermionic bubbles describe the decay of the pion in p — h and 
A — h states and take into account the conventional nuclear matter polarization effects. All these 
diagrams are responsible for the interaction of two nucleons in the particle-particle ladder channel 
with the in-medium virtual excitations. In Fig. El the first graph is the vacuum contribution. The 
second and third diagrams correspond to the p — h RPA series and the last diagrams accounts for 
the excitation of 2p — 2h states and the contribution of the 5-wave optical potential, respectively. 
The p-wave pion self energy is given by 

-l 

(10) 

where g' = 0.6 is the Landau-Migdal parameter [41] , p is the nuclear matter density and U(k, p) = 
U d {k, p) + U c (k, p) is the Lindhard function accounting for the direct and crossed contributions 
of p — h and A — h excitations with the normalization of the appendix of Ref. [42]. 
The OPEP in the momentum space takes the form 

- iV PEp(q, P) = -iW(q, p)q i q j a\a 1 2 T 1 T 2 

where we have defined 

W(q, p) = {^jf) 2 1 2 F 2 (?) Ar(0, q) (12) 
5 



n(M 



D + F 

2fn 



k 2 U(k,p) 



D + F 

2f« 



g'U{k,p) 



(11) 
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Fi gure 3: Left panel: The OPEP in the vacuum (solid curve) and at normal po — 0.16 fm 3 nuclear 
matter density (dashed curve). Right panel: OPEP at several densities. 



where D n (qo,q) is the pion propagator in the medium 

iDAk) = k*-mi-n(k, P ) + io+ (13) 

and F(q) stands for a monopole form factor A 2 /(A 2 + q 2 ) with the cut off scale A = 1 GeV, 
q is a momentum transfer and q = q/\q\. Note, that OPEP depends on the real part of the 
polarization operator only, since for q = (0, q) one has ImII(0, q, p) = 0. The well-known vacuum 
NN amplitude is recovered in Eq. (|12|) by setting IT = or in the limit p = 0. 

Our results for W(q,p) are presented in Fig. 01 (left). The standard vacuum behavior is 
shown by the solid curve. The dashed curve represents the modified OPEP at normal nuclear 
matter density were we observe an additional strong increase of the strength associated with the 
attractive excitation of p — h and A — h collective states, with p — h playing a dominant role. 
Their individual contributions are shown by dot-dot-dashed and dot-dashed curves, respectively. 

The analytic properties of the in-medium pion propagator which we use here may be verified 
by using the dispersion representation for the Green function 

iD n (k ,k) = -i- dxx-py (14) 

^ Jo K - x + ? 

In this case the OPEP of Eq. (|12|) can be written in terms of the absorptive part of the pion 
propagator only 

„ oo 
2 (D + F\ 2 N fdx 



W(q,p) = - [^j—) <T^(<?) J — ^D n (x,q) ( 1,5) 



o 



The causality requires that both equations must produce the same result. Indeed as one can see 
in Fig. |21 the curves calculated with dispersive and absorptive parts of the pion propagator are 
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Figure 4: The OPEP with short range correlations (solid curves) as a function of the nuclear matter 
density. The left and right panels are the central and tensor parts, respectively. 



practically indistinguishable. We would like to note that a strong modification of the OPEP at 
finite baryonic density observed here is not new and was predicted long time ago by Migdal [41]. 
There it was also shown that in-medium modified OPEP helps to explain the unnatural parity 
states in finite nuclei, for instance, the shift of 0~ state in closed shell nuclei. 

In the right panel of Fig. El we show our combined plot for a few densities po,po/2 and po/4- 
Here, we also show the sensitivity of our results to the value of the Landau-Migdal parameter 
g' . We find that the increase of g' from 0.6 to 0.7 makes the OPEP softer. This fact suggests 
that the proper treatment of the NN short range correlations is important for understanding the 
in-medium properties of the OPEP. 

At this point we would like to mention that in a realistic calculation one will have to add 
strong repulsive forces at short distances. This can be done in a straightforward way using any of 
many body schemes discussed in the introduction. The correlations of this part of the interaction 
would effectively modulate the 7r exchange interaction, introducing the correlation parameter 
g' [43]. The denominator in Eq. (|13|) takes into account this effect between p-wave bubbles in the 
diagrams of Fig. but not between the external nucleon and the contiguous bubble. To account 
for this we make the separation between the longitudinal and transverse parts of the pion effective 
interaction [44,45] 

{Wf m \t-^-% + io+ ^ Wq) ** + V ' i<m ' - **> (16) 

where qi is the Cartesian component of the unit vector q = q/\q\ and 

'D + F 



Q 2 



(?0 — q 2 — m\ + z0 H 

> + i 

2U 



g' 



F(qY (17) 



Vt(q) = (^-) 9'F(qf (18) 
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When we perform the sum of diagrams in Fig. [21 then we get 

Hi) ~- , v t (?) 



VopEp(q,p) 



l-W(g,p)V,(g) 



l-W(g,p)V t (<z) 



W (q, p)ct 1 (t 2 t 1 t 2 + W 2 (g, p){a- 1 qa 2 q - -o- 1 ct 2 )t 1 t 2 



(19) 



Wo(<?,p) = - 
^ 2 (<Z,p) 



+ 



2V t (g) 



l-W(g,p)Vi(g) l-W(g,p)V t (g) 
Vj(g) V t (?) 



(20) 



(21) 



l-W(g,p)VKg) l-W(g 5 p)H( g ) 

where we have explicitely separated the central and the tensor parts of the interaction. In Fig. 0] 
we show the results for the central Wq (left panel) and tensor W 2 (right panel) parts (omitting 
the spin-isospin operators) as a function of the baryonic density. Note that in the limit p — * 0, 
W 2 = W where W is given by Eq. El 



4 Two pions in the medium 



In this section we turn to the dynamics of the two pion system in the nuclear medium. Our main 
interest here is the propagation of two 5-wave pion pairs. As it is well known in the 5-wave 
scattering the use of the proper unitarization schemes lead to the generation of the er-meson. 
Later on we will use this result for the in-medium NN interaction mediated by exchange of two 
correlated pions in the scalar-isoscalar cr-meson channel. 

For the ir a TTb — ► n c iTd scattering process, defined by the Cartesian isospin indices a,..., the use 
of the standard xPT procedure in expanding the C 2 of Eq. © to order 0(tv a ) results in the tree 
level contact interaction 



iV, 



ab—>cd 



SabS c dM s ) + 6ac6bdA(t) + 8 a d5bcA( u ) 



where 



(22) 
(23) 



A(s) = j- 2 (s-Ml- 1 - Y, K)+0(q% 

n i=a,b,c,d 

and Aj = kf — M% is the off-shell part of the invariant -kit amplitude. At this order of the pion 
field expansion the isoscalar 5-wave mr partial amplitude (L = 0) is obtained from the standard 
decomposition 



V, 



LJ=0 



1 1 



2(^2> 



dcos6 P L (cos9) V^=°{9) 



(24) 



-i 



X + 

y a y 



Figure 5: The unitary series representing irir scattering amplitude. 
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Figure 6: Rcnormalization of the in- medium tttt amplitude including p-hole and A-hole excitations. 



where Pl(cos9) are the Legendre polynomials and (V2) a accounts for the statistical factor oc- 
curring in states with identical particles: a = 2 for tttt — * tttt in the unitary normalization of the 
states [31]. The tree level scalar-isoscalar tttt scattering amplitude is 



V, 



L=I=Q 



1 
72 

J 7T 



Ml 
2 



1 



3 ^ 



A, 



(25) 



In Eq. (|25|) the off shell part depends on choice of U and is equal to zero for on mass shell pions. 

Following Ref. [31] and using the Bethe-Salpeter equation we unitarize the S-wave tttt scat- 
tering amplitude (see Fig. [5J) 



V, 



L=I=0 



-,L=I=0 



Si i > 



Mi 



(26) 



where G n7T (s) is a scalar two-pion loop function 



G W7T (s) 



(IttY [(P - k) 2 - M2 + i0+](k 2 -M 2 + »0+) 



(27) 



where = {Pq,P) and s = P 2 . The G^ls) function is analytic with a cut along the positive 
real axis starting at the tttt threshold. Note that, Eq. I)26|) contains a pole in the second Riemann 
sheet corresponding to the cr-meson with mass and width M — iT/2 ~ 450 — i221 MeV. 

In the nuclear medium the S-wave tttt scattering amplitude and therefore the cr-meson get 
renormalized and explicit calculations were done in Refs. [26,46]. The diagrams at one loop 
level are shown in Fig. HO For instance, the amplitudes corresponding to the insertion of fermion 
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bubbles in the upper meson line are given by 



vr(s) 



D + F 



2f: 



i 



9/4 V 2/ w 



k 2 U(k)D 7T (k)D 7T (P - k) 
k 2 U{k)D n (P -k) 



i2 f D + F 



2U 
D + F 



d 4 k 

d 4 k 



Kn(s) + 3^^ At 



(2tt)< 



(28) 
(29) 
(30) 



As was shown in Refs. [26,46], in the center of mass frame of the two pions the off shell part (A 
terms) in V™ p cancels exactly the Vjf and Vjfj terms. And one is left only with the diagrams of 
type (I) but with the on shell tttt amplitude 



V up (s) 



D + F 

2fn 



Von(s) 



d 4 k 

k 2 U(k)Dl(k)D n (P-k) 



(27T) 



(31) 



It is straightforward to iterate the p — h and A — h excitations in Fig. H3 (I) and the loop function, 
T(s), is given by 

T(s) = V 2 n {s) ■ G n 



(s) 



where G„(s) is in-medium modified scalar loop integral 

d A k ~ 



D^{k)D n (P-k) 



(32) 



(33) 



Using the spectral representation for the in-medium pion propagators, Eq. ()13|) . we get for the 
loop function 



G WW (P) 



7T^ 



dk 

(27r)' 



dx 



ImD n (x, k) 



dy lmD n (y — x,P — k) 

x V 

(P + y)(P -y + iO+) 



(34) 



We refer to Ref. [47] where different aspects of 5-wave tttt scattering in the nuclear medium 
are discussed and also the behavior of the a-meson mass and width at finite baryonic density 
is addressed. But here, we would like to illustrate the impact of the nuclear medium on the 
tttt system. For that consider the imaginary part of the loop function in the tttt center of mass 
frame P M = (Po, 0) with Pq = ^/s. This situation is relevant for the in-medium tttt scattering and 
contains the proper information about the dynamics of the pole position of a at finite density. 
Our results for the imaginary part of the scalar loop function for several densities are shown in 
Fig. deleft panel). The solid curve correspond to the vacuum loop function which can be obtained 
from Eq. ()34|) by substituting the imaginary part of in-medium pion propagators by their vacuum 
expressions 



1 / dm 2 
ImG„(^)=--Jl * 

167T V S 



(35) 
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Figure 7: The imaginary part of the in-medium scalar loop function in the tttt CM frame (left 
panel) with the normal dressing of the pion propagators including p — h, A — h and 2p — 2h excitations. 
The right panel shows the loop function in the space like region P M = (0, q). 



which is proportional to the two-body phase space of two particles (pions). As one can see in 
Fig. (left) the effect of the medium is remarkable due to the increase of available for pions phase 
space because of additional pion decay branches like p — h, A — h and 2p — 2h. 

The kinematics relevant for the NN force, where two nucleons interact by exchange of mesons 
is defined by the moving reference frame where P^ = (0, q). In this case the tttt interaction in the 
medium has to be generalized from the results in [26,46] since there Pq ^ but q = 0. We repeat 
all the steps that led to cancellations in [26,46] and find again the same cancellations as before 
but with some remnant q 2 dependent terms vanishing in the limit q 2 = 0. To evaluate these 
terms we simplify the calculation assuming q relatively small (this is fine for momenta below the 
Fermi momentum). Concretely, we assume 

1) q 2 l^tnax ^ 1; 

where k max is the cut off in the three momentum (of the order of 1 GeV in [31]) that one uses to 
regularize the G nn function. We also assume 

2) k 2 D w (k) ~ -1 

which implies 3> m 2 as it is the case. And 3) we expand D 7r (k) in terms of D n (P — k) 

and vice- versa to relate different terms. After all this is done we find that the corrections can be 
taken into account by means of the change in the Lindhard function 



U(k) ^U{k) 



1 + 



2 



+ 



77^ 

2 



(36) 



The expression in brackets in Eq. ()36|) is ~ 1 for very small q and also ~ 1 for q 2 3> m 2 /2. Hence 
with very good approximation we can take the bracket equal to unity and thus there are no other 
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Figure 8: Diagrams representing the correlated two pion exchange with NN, NA and AA intermediate 
states. 



corrections to be done to the result of [26,46] except the obvious one of changing s — > —q 2 - For 
these value of s the G nir is only real, contrary to the case studied in [26,46]. Results for G 7T7T (q 2 ) 
can be seen in Fig. [7| (right panel) for different nuclear densities po, po/2 and po/4. We can see 
that the corrections are sizable particularly at small values of |q|. 

5 In-medium renormalization of the correlated two 
pion (cr-meson) exchange 

The correlated two pion exchange (CrTPE) in the scalar-isoscalar channel, the equivalent to a 
a exchange in meson exchange models [21], or correlated two pion exchange in the dispersion 
relations in [48] was studied in [29, 49] within the context of chiral Lagrangians. One starts 
from the diagrams of Fig. |H1 where the tttt — > irir scattering shows the off shell ambiguities. To 
avoid these ambiguities with isoscalar exchange it was stated in Ref. [49] that one must include 
the subset of diagrams of Fig. to find cancellations of the off shell tttt isoscalar amplitude. 
This statement was rigorously verified in Ref. [29], where it was shown that the consideration of 
these subset of chiral diagrams, Fig. including the contact 3ttNN interactions, results in the 
cancellation of the off-shell part of the tttt amplitude, and the on-shell part of the tttt amplitude 
can be factorized out from the loop integrals. One step forward was given in [29], where iteration 
of the tttt interaction, through the Bethe-Salpeter equation, was done by means of which a simple 
analytical expression was obtained for the correlated two pion exchange in the scalar-isoscalar 
channel 



V a {t) = 6V 2 (t) 



-fUt-^fy 1 -cut) 



(37) 



where t = —q 2 in the NN cm. frame. The vertex function V(t) = V]\r(t) + V&(t) for the triangle 
loop with two mesons and one baryon propagator, including N and A intermediate states, is 
evaluated in [29] using a cut off in \k\ of about 1 GeV. Note that, the bracket in Eq. (|37|) contains 
a pole in the s-channel corresponding to the u-meson. This restores the relation to the u-meson 
exchange, which now enters the formalism as a dynamical resonance in the tttt system [29] (see 
also related discussions in Refs. [50,51]). 

The diagrams responsible for the renormalization of CrTPE in the nuclear medium are shown 
in Fig. I1UI There, in analogy to what was done for the tttt interaction we include tt selfenergy 
corrections as well as vertex corrections. We find that the cancellation of the off shell part of the 
vacuum u-exchange discussed in [29] is also exact at finite baryonic density but for zero momentum 
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Figure 9: The set of diagrams which cancel the off- mass shell part of the scalar-isoscalar correlated 
two-pion exchange. 



transfer only (see Appendix B). The results of the derivation for g / can be summarized as 
follows: 

1) The expressions for the triangle vertex functions Vn and V A of Ref. [29] are obtained in 
the same way replacing the two free pion propagators by the renormalized ones. 

2) The Lindhard function entering the pion selfenergy is changed to 



U(k) ^U{k) ■ 



1 



1 



6 q 2 + tt4/2 



(38) 



to account for the corrections obtained at q ^ 0. These corrections are negligible for small \q\ 
and for \q\ > m n of the order of 15%, hence moderate in all cases. 

3) The final expression for the potential V a is given by Eq. (|37|) . which accounts for the 
7T7T rescattering, by substituting G nn (t) by G^{t) of Eq. and taking the expression for the 
in-medium vertex function V(t) = V/v(i) + ^a(^) where 



V N (t) 



d?k M N / 

1 k 2 + kq 



2kn 

J (2vr) 3 E(k) 



dx 



Im D^{x, k) 



dy 



x + y + E(k) - M N 



(x + y)(x + E{k) - M N )(y + E(k) - M N ) 



lmDn{y,k + q) 



(39) 



VA(t) 



2k a 4 f d 3 k M A 



k + kq 



dx 



n 2 9 J (2tt) 3 E A (k) 



oo 

x + y + E A (k)-M N 



Im D 7T (x, k) 



x / dy- 



(x + y)(x + E A (k) - M N )(y + E A (k) - M N ) L 

o 

The coupling constants K n are defined by 



JmD„(y,k + q) (40) 



K N 



D + F 



K A 



3 D + F 
V2 2U 



(41) 



There is still one more correction to be done in order to account for short range correlations. 
So far we have replaced the free pion propagator by the renormalized one of Eq. (|13|) . However, 



13 



P, Pj-k P, + q 


Pi 


Pj-k Pj + q 




Pj-k P, + q 




Pj-k Pj + q 


^ * ► f * ^ 

>< + 

/ / \ l-q 
/ \ 


— 


X, ' o"' 4 

:< + 

/ / \ l-q 
* \ 




A: \ -^+9 

V' + 

,p <!• 


— 


— * ^ j+ — w 

k \ /k+q 


P 2 P 2 ~l P 2 -q 

a) 


*2 


P 2 -l P 2 -q 

b) 


^2 


P 2 -l P 2 -q 

c) 


Pi 


P 2 -l P 2 -q 

d) 


Pj P } -k Pj + q 


Pj 


P]-k p l + q 


Pl 


P,-k P] + q 


Pi 


P 7 - P 2 + ? 


> r * T~* — 

k \ It' k+ v 

> — > — — 


— >- 


~r * r - * — 

"& + 

> v ' > 


—>- 


* — 


—>- 


* * — 

' ^ , 


P 2 P 2 ~ I P 2~ 1 

e) 

p 2 p 2 + ? 


— >- 

Pl 


P 2 -/ P 2 -9 
f) 

P ; P.-fc P i + <? 


— *- 

Pl 

Pl 


P 2"' P 2 "? 
g) 

P,-A: P i + 9 


— >■ 

Pi 

Pi 


P 2 -l P 2 -q 

h) 

Pj-k Pj+q 


I / \l~q + 




AO 7 *** 

l / \ l-q + 


I 


— > 

k \ /k+q + 




— * — 

k \ / k+q 

^ — — — 


p 2 p 2 -i P 2 -q 
i) 


Pl 


P l~ l Pi-q 
j) 


Pl 


P 2~l P 2 ~q 

k) 




p i p r l P 2 ~q 

1) 



Figure 10: The in- medium diagrams involving the p-wave pion dressing and vertex corrections. The 
diagrams with the intermediate A states are not shown. 



since the interaction between the p — h bubble and the external nucleons is affected by correlations 
one should take Vi (see Eq. 1)17)0 instead of [(D + F) 2 /2f. K ] 2 D 7r . Thus we would get the series 

D (k) 

D v (k) + D v (k)U(k)Vi(k) + D v (k)U(k)Vi(k)U(k)Vi(k) + ■ ■ ■ - T> ' 



l-U(k)Vi(k) 
D n (k) 



(42) 



2/. 



g'F z (k)U(k) 



where D w and D n are the free and dressed pion propagators, respectively. Hence, the expressions 
for Vn and Va get modified by including inside the J d 3 k integral the factors 

1 (43) 



^ l (k)U{k) l-(-^_J g>F\k + q )U{k + q) 

Our results for cr-meson exchange in the momentum space are shown if Fig. ^Jfor both cases, 
no short range correlations (left) and with short range correlations (right). 
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6 Uncorrelated two pion exchange 

In this section we consider another sort of intermediate distance contributions to the NN force 
generated by the uncorrelated two pion exchange. The material presented here for the vacuum 
NN scattering is standard and we merely generalize it to the nuclear medium. 

In the perturbative expansion of the NN force we must to take into account the planar and 
crossed box diagrams shown in Figs. H2l and respectively. We will discuss the scalar-isoscalar 
part of this contributions only. The contribution of the isovector exchange is small and can be 
found, for instance, in Ref. [52] 

In vacuum, the expression for the planar box diagrams, Fig. 1121 including the nucleon pole 
and NA and AA intermediate states reads 



-*V, 



(P) 



NN 



d 3 k 
(2^)3 



E« • (k + q)] [SJW • k] 1 [Eg) • (k + q)] [Eg? • k] ] {ni p 2(k, q) } 

(44) 

In Eq. (|44p the sum over n, m = N, A is assumed and the spin transition operators are E^r = <r, 
Sa = S. The coupling constants n n are defined in Eq. lHJ) . The isospin factors are given by 



r(P) 



for which we find 



T (P) 

NN 



3-2tit 



1^2, 



l NA 



2 + -TiT 2 , 



r(P) 

'aa 



Ti To 

3 9 



(45) 



The function ?lS in Eq. Q44|) contains the integration over the time-like component of the four 
vector k. 



dk^ 



D n (k + q) D„(k) 5 n (Pi - fc) 5 m (P 2 + k) 



(46) 
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Figure 12: The planar box diagrams involving the nucleon and A intermediate states. 



where D w and Si is the pion and nonrelativistic baryon propagators, respectively. Si is given by 

Mi i 



iSi(P) 

a/P 2 + M 2 P - X /P 2 + M 2 + i0+ 
This integration can be carried out explicitely 



(47) 



i 

2 L 



uj\ + uo 2 2 + 3^ 2 + [E n - E] [E m - E] ) (E n + E m - 2E) 

+ [uji+ w 2 ) {2UJXUJ2 + [E n + E m - 2E} 2 ^ 
1 f 1 



ljil<J2(wi + W 2 ) i -E n + # m - IE - iO H 
X E n \u x + E n - E) [lu 2 + E n -E] 

M m 

E m [u\ + -E m — E] [u>2 + S m — E] 
Recall that g = (0, q), P x = (E, P), P 2 = (E, -P) in NN cm. frame and in Eq. (gHJ) 



(48) 



^ - X /fc 2 + M2, W2 = ^(k + qf + Ml E n = y/(P-k) 2 +Ml E m = y/(P - fc) 2 + M m , 

(49) 

are the on-shell energies of intermediate pions and baryons with cm. energy of the initial nucleons 
E = \/P 2 + m. 

In the nuclear matter 1Z — > TZ, and the expression for 1Z is obtained by using the dispersion 
representation for the in-medium pion propagator 



oo oo 



ftnm( fc >9»P) = ^2 J dx 2 j dy 2 TZ ( n P 2(x,y)^D 7T (x,k, p)ImD n (y,k + q,p) 



(50) 
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Figure 13: The crossed box diagrams with nucleon and A intermediate states. 



where 1Zrfm{x,y) is given by Eq. with substitution oj\ — > x and u; 2 — ► y. 
For the crossed box diagrams, Fig. I13( we have 



_ jV (C) _ V K K 

NN — /, "fiAm 
n,m 



d 3 k 



(2tt) 3 

The isospin factors are the following 

r(C) 



£ £> • (fe + g)] [e« • fc] 1 [e$ • fe] [e£> • (fc + g)] 1 (fc, g) } 

(51) 



/^ = 3 + 2TiT 2) I^l = 2-^ Tl T 2 , J^=g + gTiT2 

and the vacuum expression for 1Zml{k, q) is given by 

00 

7d(fc,g)= y ^^(fc + g)^^)^^- fc)S m (P 2 -A: 

—00 

The integration can be done analytically and our result reads 



r(C) 



4 2 



(52) 



(53) 



w 2 



uf + uj + LO\Ul2 + (Wi + W 2 )(£„ +^m- 2£) + (£ n - E)(E m - E) 
1 



U)\U)2{U\ + OJ2) 



Mr, 



E n [u 1 + E n - E\[u2 + E n - E\ 



E m [ui + E m — E] [u 2 + E m — E] 
here W\ , oj^-, E and E m are defined in Eq. (|49() and 

£ n = y/(P + k + qf + M r ] 



(54) 



(55) 
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The corresponding expression for the in-medium crossed box diagrams takes the form 

oo oo 

^nm( fc 1 g>/') = ^2 J dx2 j d y 2 K^(x,y)Imbn(x,k,p)lmDi r (y,k + q,p) (56) 
o o 

For the generic case of non- vanishing initial momenta the analytical structure of Eq. ()48l) is 
driven by the term in figure brackets 

' (57) 



E n + E m -2E- iO 



Here we have to pay a special attention to the case where two intermediate NN states n = m = N 
appear, because in time ordered perturbation theory these diagrams are generated by iterations 
of the OPEP in a Lippmann-Schwinger equation (LSE). Considering the NN intermediate state 
only, the iterated TPE can be easily identified and comes from the nucleon pole in Eq. (|46|) (in 



the lower half of the complex plane) corresponding to ko = En — E — i0 + with E = y P + M^- 

1 1 



/ \ 2 

^{P)N-pole (i , , ^ _ — if Mn\ i 



{uJuUJ2) -Y{e^) E-E N + iO+ 



(E-E n )*-lj* {E-E N f-ul 



(58) 



Expanding Eq. ()58j) in powers of 1/M 



v 



* ->^ p l_ (p M _ N kr+i0+ + 0(1/M N) (59) 

and inserting the leading order result in Eq. (|44[) one can get the second order term in the 
non-relativistic Lippmann-Schwinger equation. 
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It is instructive to derive the contributions of the two remaining poles (in lower half-plane 
from the pion propagators 

E N ? 



n 



(P)-k— poles 
NN 



(k,q) 



Uj\ + U>2 + 0J\Ld2 



After the expansion of this result in powers of 1/Mjv we get 



eg 

E N f)(ul 



(E - E N f 



Ml 

E N 



K 



(P)tt— poles 
NN 



(k, q) ~ 



2 uf^i^i + 0J2) 



+ 0(1/M N ) 



(60) 



(61) 



In this limit Eq. Q61JI cancels exactly the corresponding crossed box diagram with two interme- 
diate nucleons in the isoscalar channel. Indeed, considering the crossed box diagram with two 

(C) 

intermediate nucleons, the leading term of TZ N ^(k, q) in the 1/Mn expansion is given by 



K {C) (k a) ~ i <d±4 



UJ1UJ2 



2 LofiO^iuJl + LO2) 



+ 0(1/M N ) 



(62) 



Note that for the isovector tttt exchange because of the different sign of the T\T2 term in I 



AC) 

NN 



(P) 
NN ' 



in Eqs. ()45l) and <\52\i these two contributions would add. The cancellation discussed above 
hold also in the nuclear medium. 

The result of the vacuum scalar-isoscalar iViV force generated by the planar and crossed box 
diagrams, with the nucleon pole diagrams excluded, is shown in the left panel of Fig. El by the 
solid curve. It is in agreement with Ref. [52] where it was shown that the consistent use of the cut 
off regularization in both the correlated two pion exchange, and uncorrelated two pion exchange, 
together with the contribution of a repulsive w-exchange lead to a scalar-isoscalar potential in 
good agreement with the Argonne [53] potential in the whole range of relevant distances. The 
corresponding results for the nuclear matter are shown in Fig. ^] (left) for three densities po, 
Po/ 2 and po/4- Qualitatively the behavior is similar to the vacuum case but with a strong 
enhancement toward the small momentum transfer. As we have seen this feature is generic in 
present calculations. Uncorrelated tttt exchange becomes extremely attractive at intermediate 
distances. 

The spin sum over intermediate baryon states of the ~ 53 • k operator in Eqs. (|44|) and (|51|) 
gives in the scalar channel 



P[(k + q)-kY 



(63) 



where = 1 (NN), (3 = 2/3 (AT A) and /? 
the short range correlations 



4/9 (AA). Now we again wish to take into account 



D + F 

2fn 



F 2 (k)F 2 (k + q)D„(k + q)D n (k) [(k + q) ■ k]< 
=^ Wi{k + q)Wi(k) Uk + q) ■ k 



+ 



m(k + q)Wt(k) + Wi{k)W t {k + q) 



(k + q)-k 



+W t (k + q)W t (k) I 1+ (k + q)-k 



(64) 
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Figure 15: The modification of the NN interaction at normal nuclear matter density in momentum 
space (left panel) and in configuration space (right panel). Shown are the central part of the OPEP (solid 
curve), tensor part of the OPEP (dashed curve), correlated a- meson exchange (dot-dashed curve) and 
uncorrelated irir exchange (dot-dot-dashed curve). 



where 

*<*>-dlb) <65) 

where U{k) is the Lindhard function and Vi and Vt are given by Eq. (fTTj) and (|18|1. respectively. 



7 Results and discussions 

By comparing the results in Fig. and Fig. 0] we observe that the effect of correlations has been 
essential and reduces drastically the medium effects found in Fig. |3] without corrections. We 
observe that the medium corrections weaken the strength of the central part of the OPEP. On 
the other hand the effect of the medium corrections on the tensor part are more moderate. 

The a exchange presents similar features as we can see in Fig. 1111 There we see that the 
medium effects in the absence of short range correlations are rather large and increase the strength 
of the interaction in about a factor two at p = po- However, as soon as the correlations are taken 
into account the medium modifications are reduced drastically and at p = po just reduce the 
strength in about 25 percent. 

The medium effects in the uncorrelated two pion exchange are also relatively small, with 
respect to the size of the interaction in the vacuum, see Fig. 14 (right). Once again the consider- 
ation of the correlations has been essential and reducing the size of the medium effects. However, 
the fact that the strength of this interaction is bigger than that of the correlated two pion ex- 
change makes the absolute correction relevant and of the same size as the corrections discussed 
before. 

Since we have taken also the effect of correlations in the free part of the interaction and any 
realistic calculation of the binding energy of matter will explicitly account for these correlations, 
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the use of our modified in-medium potentials would lead to double-counting if any of these 
methods is used. For this purpose the relevant results from this work should be the differences 
between the in-medium potential and the free one. In this case we start already with one bubble 
and the correlations account for the repulsion between the external leg and the one in the bubble. 
The two external legs still have to be correlated and this will be done with the use, for instance, 
of the ordinary Bethe Golstone equation. 

After this discussion of the medium effects in the different terms, we show the differences in 
momentum space between the medium and free parts of the different terms in Fig. ^3 (left). The 
results are shown at p = po and we see that these corrections are moderate, but they could have 
a relevant role in the binding of matter. In order to have a qualitative idea of the relevance of 
these corrections to the potentials we rewrite them in coordinate space and show these results 
in Fig. El (right). The effects seem sizeable at short distances, but this will be irrelevant in any 
realistic calculation of binding energies since the consideration of the short range correlations 
between the external legs will make this interaction inoperative. More interesting is the strength 
of the corrections around 1 fm, and there we see that all them are relatively small, of the order 
of 20 MeV or less, the biggest one being the central part of the OPE. However, given the size 
of the empirical scalar isoscalar attraction [53] which is of the order of 20 MeV at intermediate 
densities, the corrections found here are not negligible. It would be thus interesting to see the 
effects of the results obtained here in observables like the binding of nuclear matter and other 
properties, which we hope to stimulate with the results obtained here. 

8 Conclusions 

In this paper we studied the modification of the one pion exchange, as well as two pion exchange 
potential inside a nuclear medium. For this purpose we separated the two pion exchange into 
an uncorrelated two pion exchange and the correlated two pion exchange. We study both in the 
scalar isoscalar channel, which is by large the most important one generated by this interaction. 
The correlated two pion exchange gave rise to the equivalent of the a exchange in other models 
and we studied the medium modification to it. On the other hand, for the uncorrelated two 
pion exchange we have followed a traditional approach in which only terms with at least one 
intermediate A state are considered. 

One of the important findings here was the effect played by the NN short range correlations 
which drastically moderated the medium corrections to the potential. In the absence of these, the 
corrections where unrealistically large. Yet, even if relatively moderate, the medium corrections 
found in this paper are sizeable enough to have relevant repercussions in the binding and other 
properties of nuclear matter and we would like to encourage calculations in this direction. 
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A Feynman rules for vertices 



The set of Feynman diagrams shown in Fig. ED appear in the construction of the irN — > ttttN 
transition amplitude. The corresponding 3ttNN vertex functions are given by 



\ 7t 71 V' , \ 71 TZ , \ 71 7C °/ 

2 \ /<ii 2 \ yV 2 2 *\ 




Figure 16: Contact terms appearing in the construction of the irN — > irirN transition amplitude. 
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B Correlated exchange in the medium 

In this Appendix we demonstrate the cancellation of the off shell part of the correlated two pion 
exchange in the scalar isoscalar channel. For simplicity we consider the nucleon intermediate 
states only. The generalization to A^A and AA is straightforward. The diagrams are shown in 
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Fig. EI] and corresponding amplitudes are given by 
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where Va^tt^) is the triangle loop integral 
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and the on-mass-shell tttt scattering amplitude is 
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In the present case g = (0, q) and for q — ► the amplitude cancels exactly Tj. The same 
>n is found fo 

- iTP - iT P 



cancelation is found for and Tj. The sum of T e and is given by 
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The sum of Tf and takes the form 
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From this the sum of four diagrams is given by 
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Alternatively, the spin flip parts of Eqs. (|§U|) and (|§Tj) can be canceled if we change in Eq. (|§Tj) 
k — * —(k + q) and using again the fact that q = (0, q) and the properties of the Lindhard function 
U(p Q ,p) =U{po, -P] we get 
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and 
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Finally 
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As one can see in the limit q — > we get an exact cancellation of the off mass shell part of the 
-kit amplitude and only T° n and T£ n are left. 
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